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Abstract 

A separable potential formalism is used to describe the tttt and KK 
interactions in the I G (J FC ) = + (0 ++ ) states in the energy range from 
the tttt threshold up to 1.4 GeV. Introduction of relativistic propaga- 
tors into a system of Lippmann-Schwinger equations leads to a very 
good description of the data (x 2 = 0.93 per one degree of freedom). 
Three poles are found in this energy region: /o(500) (M = 506 ± 10 
MeV, r = 494 ±5 MeV), / (975) (M = 973±2 MeV, T = 29±2 
MeV) and / (1400) (M = 1430±5 MeV, V = 145 ± 25 MeV). The 
/o(975) state can be interpreted as a KK bound state. The /o(500) 
state may be associated with the often postulated very broad scalar 
resonance under the KK threshold (sometimes called a or e meson). 
The scattering lengths in the tttt and KK channels have also been 
obtained. The relativistic approach provides qualitatively new results 
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(for example the appearance of the /o(500) ) in comparison with pre- 
viously used nonrelativistic approach. Interactions in both channels 
are attractive and have short range form factors. 

PACS number(s): 14.40.Cs, 13.75.Lb, 12.40.Qq, 11.80.Gw 



I Introduction 

Scalar meson spectroscopy is still far from being well understood [1-4]. 
Since it was difficult to explain known properties of the l G (J PC ) = + (0 ++ ) 
states using a standard qq picture, other models treating the four-quark states 
||, H| or meson-meson molecules [7-11] have been invented. There is also a 
continuous search ]T2l IT3J for scalar gluonium states which can be mixed with 



the quark states of the same quantum numbers. The nature of the /o(975) 
and ao(980) mesons is very controversial [14-17]. Closeness of their almost 
degenerated masses to the KK threshold energy constitutes an argument 
towards their interpretation as virtual bound KK states which are unstable 
due to the open tttt channel [10, 18-21]. 



In this article we extend the coupled channel formalism of Refs [y| [19 
including some relativistic effects in both the nn and KK decay channels 
of the scalar mesons. We use the relativistic propagators and the separable 
potentials in the Lippmann-Schwinger formalism. Such an approach has 
already been applied in the analysis of the pion-nucleon amplitudes p2[ or 



in the studies of the A (1405) resonance structure [2J]. Our aim is to 
describe quantitatively the isoscalar s-wave tttt and the KK scattering data 
in a wide energy range starting from the inr threshold up to 1.4 GeV. The 
parameters of the meson-meson interactions are fitted to the data and the 
S'-matrix structure is analysed in order to extract information about the 
scalar resonances in that energy region. This is a new step in comparison 



with [pl| where some pole structure of the S'-matrix has been postulated 
and the \ 2 tests have not been done. We observe that the relativistic effects 
are important not only in the nir channel but also in the KK channel even 
near the threshold. This fact may have important implications, for example 
the /o(975) meson, interpreted as a KK molecule, can still have a relatively 
small radius (less than 1 fm). Consequently, the predictions of the radiative 
(j) decay into the KK 7 system via <p — > fo (975)7 have to be substantially 
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influenced Pq| . The importance of the relativistic effects on the two photon 
decays of the J PC = ++ and 2 ++ states has also been recently stressed [^6| . 

This paper is organized as follows. In Sect. II we define our formalism. 
Sect. Ill is devoted to an analysis of the single tttt channel essentially up to 
the energy of about 700 MeV. The properties of the KK channel interaction 
in the decoupled case are discussed in Sect. IV. Both the nonrelativistic and 
relativistic cases are treated. In Sect. V we analyse the interactions in the 
coupled 7T7T and KK channels. A comparison with the data is performed 
and the S-matrix structure is studied. In Sect. VI we summarize the main 
results and outline some perspectives. 



II The Formalism 



We describe the tttt and KK isoscalar s-wave interactions in the frame- 

The scat- 



work of the coupled channel Lippmann-Schwinger equations |22 
tering amplitude T satisfies the following equation in the momentum space 



< P I T | q > 



< P I V | q > + 

d 3 s 

- — — < p I V I s >< s I G I s >< s I T I q >, (1) 



where V, G, T are 2x2 matrices (label 1 denotes the KK channel and label 2 
- the hit channel), V is the interaction matrix and G is the diagonal matrix 
of propagators written in the center of mass system: 

< s \ Gij \ s >= Gi(s) 5 id (i,j = l,2), (2) 

and Gi\s) = E - 2E l (s) + ie, e — ► 0(+). I n Eq. @ s is the relative 

momentum, E is the total energy and E^s) = \J s 2 + mf are the relativistic 
energies; m!=495.69 MeV, m 2 = 137.27 MeV are the average kaon and pion 
masses. 

A very convenient parameterization of the interaction matrix is a sep- 
arable potential form. In the KK channel we use the simplest rank-one 
potential in the momentum space 

< p | V n | q >= An 9i(p)gi(q) (3) 
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and in the tttt channel a rank-two potential 



< p | V 22 | q >= A 22 02(p)02(g) + X 33 g 3 (p)g 3 (q). (4) 

The transition potential matrix elements read 

< P I V12 I q >=< q I V21 I p >= X 12 gi(p)g 2 (q) + X 13 gi{p)gz{q)- (5) 

In Eqs. © _ (H) Xik (i, k — 1,2, 3) are the coupling constants and gi are the 
form factors which we have chosen in the Yamaguchi form [27| : 



where 

m =i mK if i=1 m 

7711 \ m n if i=2,3 [t) 

and /3j are the form factor range parameters. The potentials (§^J) are nonlo- 
cal. After this choice the potential matrix has eight parameters (five coupling 
constants and three range parameters) which should be fixed by a comparison 
of the theoretical amplitudes with experimental data. 

The separable form of the interactions enables us to write the following 
scattering matrix elements: 

<p|T n |q> = 01 (p) tu & (g), (8) 
< P I T 22 I q > = 02 (p) t 22 92 (q) + 02 (p) t 2 3 03 (q) 

+ 03(^)^32 02(g) +03 (p)*33 03(g), (9) 

<p|T 12 |q> = < q I T 2 i I p >= 0i(p)ti2 02(g) +0i(p)ti3 03(g), (10) 

where t^^ are energy dependent reduced amplitudes (/c, / = 1, 2, 3). In Eq. (@) 
Tn denotes the iTii' elastic scattering amplitude, T22 in Eq. (^]) is the nir 
elastic scattering amplitude while Tyi and T21 are the transition KK — > 7r7r 
and 7T7T — >■ i^K amplitudes. The system (0) of the coupled integral equations 
satisfied by the T k[ elements leads to a set of algebraic equations for the 
reduced amplitudes t^i written in the 3x3 matrix form: 

t = X + XIt. (11) 
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In Eq. ( pTTJ ) A is the symmetric 3x3 matrix of the coupling constants 





/ A n 


Al2 


A13 


A = 


Al2 


A22 







V A13 





A33 



(12) 



and I is the following symmetric matrix 



hi 











-^22 


-^23 





/23 


-^33 



(13) 



consisting of the integrals 



hi 



dh 



■ gi(s)Gx(s)gx(s) 



and 



l kl 



d 3 i 



g k (s)G 2 (s) gi {s), (fc,Z = 2,3). 



(14) 



(15) 



In Eqs. (|5|), (|T2|) and ( |T3"D label 1 refers to the KK channel and labels 2 and 
3 to the 7T7T channel in which the interaction potential contains two terms. 
The matrices t and I are functions of the energy E (see Eq. ([2])). 
A solution of Eq. fliTD is straightforward: 



t 



[i-xi)- l x 



(16) 



and the resulting t- matrix is symmetric (£&/ = tik) since the matrices A and 
/ are symmetric. Explicit expressions for the matrix elements t^i are given 
in the Appendix. All the functions t^i are inversely proportional to the Jost 
function 

D(E) = det(l - XI). (17) 

Let us notice that the coupling constants Ay defined by ©-(§) have the 
dimension (MeV/c) 3 and the integrals (|l~4l--[l5|) - the dimension (MeV) -3 . 



For simplicity in further calculations we will use the dimensionless coupling 
constants defined as 



A 



Xij 



ij 



(18) 
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and redefined integrals 

J ij = 2(J3 i /3 j ) 3 ^I ij . (19) 

If the energy E is higher than the KK threshold mass we can re- 
late (compare JEJ) the on-shell scattering matrix elements T k[ (ki, k 2 ) to the 
^-matrix elements 

ik 1 E 1 {k 1 ) 

on = 1 7, T n {ki, k 2 ), (20) 

ik 2 E 2 (k 2 ) rr n i \ mi\ 
S 22 = 1 T 22 {k h k 2 ), (21) 

Z7T 
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S 12 = S 21 = - — Jk 1 E 1 (k 1 )k 2 E 2 (k 2 )T 12 (k 1} k 2 ), (22) 

where the KK channel and tttt channel momenta k\ and k 2 are defined by 
the energy conservation condition: 



E = 2^1 k\ + m 2 K = 2^ k\ + m\. (23) 

This relation allows one to write the functions I\\{E) and Iki(E) appearing 
in Eqs. (|b|) and (|T5| ) as functions hi(ki) and Iki{k 2 ). Similarly the Jost 
function D(E) given by ([[7]) can be expressed as a function of two related 
variables 

D(h, k 2 ) = DK^Dnikt) - C(k h k 2 ), (24) 



where 



D K {k x ) = 1-An J n {h), (25) 

D n (k 2 ) = 1 -A 22 J 22 (k 2 ) -A 33 J 33 (k 2 ) + A 22 A 33 d(k 2 ), (26) 

C(h,k 2 ) = Ai 2 Jn(A;i)[Ai2 J 22 (k 2 ) + A 13 J 23 (k 2 ) - A 12 A 33 d(k 2 )} 

+ A13 Ju(A:i)[Ai 2 J 2 3(A; 2 )+Ai3 J 33 {k 2 ) - A 13 A 22 d{k 2 )\ (27) 



and 



d(* 2 ) = J 22 (k 2 ) J 33 (k 2 ) - J 2 2 3 (A; 2 ). (28) 

In Eq. ( PB] ) Dx(ki) is a part of the Jost function corresponding to the inter- 
action in the single KK channel, Dft(k 2 ) is a similar part in the tttt channel 
and C(ki,k 2 ) comes from the interchannel coupling. 
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Using the Jost function D we can also express the diagonal S'-matrix 
elements as the ratios: 



Sn 



D(k u k 2 ) 
D(k h k 2 ) 



(29) 



(30) 



and the nondiagonal matrix element satisfies the equation 



r.2 



D{-K-k 2 ) 
D(k u k 2 ) 



(31) 



Above the KK threshold the S'-matrix can be parameterized in terms of the 
inelasticity parameter rj and the phase shifts 5^ and 6 K k'- 



If the energy E is lower than 2rriK then we can still use Eqs. (0) and (|30|). 



provided k\ is purely imaginary k\ = i\jvr? K — E 2 /A. In this case the inelas- 

ticity parameter rj = 1 since the KK channel is closed. 

Energy dependence of the experimentally measured quantities 5 nn , S K ^ 
and rj is closely related to the analytical structure of the Jost function in the 
complex planes of the k\ and k 2 momenta. 

Ill 7T7T channel interactions 

In this section we consider interactions in the single 7T7T channel i.e. with- 
out couplings to the KK channel (A12 = A13 = 0). We study the energy 
range from the rnr threshold up to about 700 MeV, where in the first step 
the influence of the coupling to the KK channel can be neglected. The 
aim of such an analysis is an estimation of the ttti potential parameters and 
investigation of the pole structure of the ttti scattering function S 22 . Two 
experimental facts can be quantitatively described in such an approach: the 
first is a positive value of the irn s-wave scattering length in the 1 = 
s-state, and the second is a systematic increase of the 1m phase shifts from 




(32) 
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degrees at the tttt threshold up to about 70 degrees at the energy of 700 MeV. 
Experimental data indicate the existence of the /o(1400) meson. Although 
the /q(1400) mass is higher than the energy limit under consideration, this 
resonance should be included in our analysis because of its large coupling to 
the 7T7T channel and a large width. 

At the beginning let us discuss when the propagator G 2 (see Eq. (0)) 

has the nonrelativistic form G 2 (s) = (k^ — s 2 ) / m 2 +ie. We start by assuming 
that in the potential V 22 in Eq. (|j) only the first term, i.e. the coupling 
constant A33 = 0, is present. Then the Jost function Djc(k 2 ) can be expressed 
as (see Eqs. (p6|) and ( |A13| ) in Appendix) 



= 1 + (T^P < 33) 

where a 2 = k 2 / (3 2 . Its phase is directly related to the phase shift 5 wn by 

D n (k 2 ) =| D n (k 2 ) I e~ iS -. (34) 
Knowing that the scattering length 

a n7T = lim 6 m /kz, (35) 

k'2 >0 

we can evaluate a low energy limit of the Jost function ([33]). As a result we 
obtain 

°- = -Ad + A*)- (36) 

The positivity of the experimental a nn value imposes the following conditions 
on the coupling constant A 22 : 

- 1< A 22 < 0. (37) 



The analytical structure of the Jost function in Eq. (|33D is very simple: it has 
a double pole at k 2 = —i(3 2 and two single roots at k 2 = i/3 2 (±y / —A 2 — 1). 
For the A 22 satisfying the inequalities (|37|) both zeros lie on the negative 
part of the imaginary k 2 -&xis. Since the matrix element S 22 is given by a 
ratio D*ji(—k 2 ) / D<ji{k 2 ) , its poles coincide with the Jost function zeros of the 
denominator except of the numerator double pole at k 2 = i(3 2 . The existence 
of the /o(1400) meson must manifest itself by the presence of two complex 
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conjugated S-matrix poles lying off the imaginary k 2 -cLX.is. It has been shown 
in [18] that the introduction of the second term in the potential V nn permits 



one to obtain such a pair of poles and, in addition, a positive value of the 
scattering length. The coupling constants A 2 2 and A 33 in had opposite 
signs. Although it was not shown there, we can prove that in the case when 
the tttt potential contains in Eq. @ two terms, when there is no coupling to 
the KK channel, and the scattering length is positive, the two other poles 
of the S 22 can lie only on the imaginary axis. 

Next we discuss the relativistic form of the propagator G 2 and investigate 
the analytical structure of the matrix element S^^ calculated for a single term 
in V n7T (A33 = 0). The situation now is quite different: the S 22 function has a 
pair of poles off the imaginary axis and simultaneously the condition a n7T > 
can be fulfilled. The expression for the scattering length in the relativistic 
case is 

2A 22 



02 



A 



22 



if P2 I m 2 j 

TT \ 777.2 "^2" 



(38) 



In order to 



where d 2 is a constant defined in Appendix by ( |A12| ) and 
get a n7T > the coupling constant A 2 2 must satisfy the following conditions: 



'5 + i (- + T* 

2 TT V 777 2 ?2 , 



< A 22 < 0. 



(39) 



Negative sign of the coupling constant A 2 2 means that the interaction is 
attractive. 

Using ([34|) we have compared our theoretical predictions for the tttt phase 
shifts with the experimental data in the energy range from the tttt threshold 
up to 700 MeV. We have used the data from the K ei decay [p8j (5 points) 
and from the reaction — ■> 7r~7r + 77 [29-31] (16 points up to 700 MeV). 
We have performed the \ 2 test using the CERN program MINUIT and as 
the output parameters we obtained A 2 2 = —0.213 and (3 2 = 1239 MeV. 
For these parameters x 2 = 1-13 per one degree of freedom. The calculated 
scattering length a n7T = 0.18 m" 1 can be compared with experimental values 
(0.28 ± 0.05) 777; 1 @, (0.24 ± 0.09) 777- 1 |2| and (0.207 ± 0.028) m' 1 [0 . 
In a combined analysis [Q of the ttA^ — * ttttN data the value a n7T = (0.197± 
0.01) m" 1 is quoted. The scattering S'-matrix function 622 has only one pair 
of the complex conjugated poles at the energy E = 491 ± 7246 MeV. It is 
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clear that it cannot be related to the /o(1400) meson which according to the 
Particle Data Group |34| has the energy 1400 — 2 (75 -j- 200) MeV. In order 
to include the / (1400) state in a further analysis we have added a second 
term to the tttt potential and the Jost function D^{k2) obtained the form 
(|26|). Then in the complex momentum plane we have found a second pair of 
the complex conjugated poles of the S22 function which can be related to the 
/o(1400) meson. In the numerical calculations of the pole positions we must 
take into account the energy range above 1 GeV and include the coupling to 
the KK channel. Therefore a full discussion concerning the /o(1400) meson 
is presented in Section V. It is worth-while to note, however, that contrary 
to the nonrelativistic case, both coupling constants A22 and A33 have received 
negative signs. 

The appearance of a new S22 pole, in comparison with the nonrelativistic 
case, at an energy about 500 MeV can be related to the existence of the so- 
called a or e meson which has been often postulated both experimentally and 
theoretically (see references in p5|). The mass of this meson varies between 
about 500 MeV and 1000 MeV, while the width even more - between 300 
MeV and 1000 MeV. Various models of the a meson structure have been 
formulated describing it as normal qq, qqqq or gluon-gluon states. Such a low 
mass scalar meson also appeared in the so-called cr-model and in other field 
theoretical models (see [^] and references cited therein) in which its mass 
had a value of about 500 MeV and its large width varied from 300 MeV to 
more than 500 MeV. The exchange of the a meson has been used in the Bonn 



model of the nucleon-nucleon interactions [57]. It has been found that the 
a meson exchange can effectively describe the contribution of two correlated 
s-wave pions in the isospin state. The a mass and width have been also 



recently discussed in the Nambu-Jona-Lasinio model |]38| . An experimental 
measurement of the a parameters is very difficult because of its large width 
[35-39]. This meson does not appear as a typical Breit-Wigner peak and in 
the partial wave analyses it can easily be interpreted as a background. The 
last time when the a meson was included in the Particle Data Group Tables 
was 1974 with M a < 700 MeV and T a > 600 MeV (and it was called e) @. 

In the isospin 2, spin tttt channel one does not expect a resonant be- 
haviour of the scattering amplitude. Using a simple rank-one potential of 
the Yamaguchi form we have obtained a very good description of the corre- 
sponding 7T7T phase shifts up to about 1 GeV. No poles of the 1=2 scattering 
amplitude have been found in the momentum range of Rek2 between and 
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4 GeV and Imk2 between and -4 GeV. This fact, which confirms the above 
expectation, is in contrast to the situation in the isospin channel in which 
we do find resonances. 

IV KK interactions 

Let us now discuss the KK interactions without coupling to the nn 
channel. We shall be especially interested in the energy region close to the 
KK threshold near 1 GeV. In this paper we do not distinguish the K°K° and 
K + K~ thresholds so we use their average energy value E t h = 2mi = 991.38 
MeV. The /o(975) meson energy is very close to this value, so it is quite 
natural to interpret this state as a bound KK state with a binding energy 
of about 16 MeV. We do not, however, postulate the existence of the bound 
state from the beginning. As we shall see in the next chapter, the parameters 
of the KK interaction will be fixed by a global fit to the tttt and KK data. 
Here we discuss different conditions for the existence of the particular KK 
structures: bound states, antibound states (sometimes called virtual states 
or resonances. For each case we shall examine effects of the relativistic 
propagation of kaons and their influence on some observables, like scattering 
length or phase shifts. A very good method to distinguish the three types of 
the KK states is to look at the positions of the S'-matrix poles corresponding 
to the Jost function zeros in the complex /ci-plane (Eq. |29|). If the zero lies on 
the positive imaginary axis (ki=iai, a>i > 0), then a bound state exists. The 
antibound state corresponds to a zero lying on the negative part of imaginary 
axis («i < 0). If there are two complex conjugated zeros in the lower half 
plane then there is a resonance in the KK channel. 

At first we analyse the most important features of the bound state. Its 
binding energy can be defined as 

E B = 2m K — 2\Jm 2 K — aj for a\ < m K . (40) 

The Jost function Dxiki) in Eq. ([25|) vanishes if 

An = l/ReiMiat)} (41) 

since Im[Ju{iai)} = 0. The coupling constant depends on the range pa- 
rameter (3\ which can be chosen regardless of the binding energy Eb- In the 
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nonrelativistic case we approximate the propagator G\ in Eq. (0) by 

G 1 (s) = ^- 2 . (42) 

Then according to Eq. ( A13| ) the dimensionless KK coupling constant is 
particularly simple 

A^-( 1 + |) 2 . (43) 

The KK force is attractive and the coupling constant is negative and smaller 
than —1. For the relativistic propagator we have numerically verified (see 
also (Eq. |A9|) ) that the /^-dependence of An is very similar to the one in 
the nonrelativistic case for the small values of Pi up to about 200 MeV. For 
larger /3i the absolute value of the "relativistic" coupling constant is smaller 
than the coupling constant in the nonrelativistic limit. If Pi — > oo, then 
nonrelativistically A n — > —1, while in the relativistic case An — > as 
—tiuik/Pi (see Eqs. (|A9|) and ([41]) ). This is the first indication that in the 
limit of large Pi the relativistic effects might be important. 

In the absence of the coupling to the tttt channel we can use similar 
relation between the KK scattering phase shift 8 K j^ and the phase of the Jost 
function Dx{ki) as in Eq. fl34]) for the tttt channel. In the limit h\ — > the 
KK scattering length is given by a K -j^ ~ ^{^kk) I ^l- The energy dependence 
of the phase shifts can be directly calculated using Eqs. fl23|), (p5l), ( |A9|) or 
(|A13| ) for a specific choice of the bound state energy at a given value of the 
range parameter f3\. This is illustrated in Fig. la for two values of /?i=150 
MeV and 2000 MeV. At the KK threshold the value of the phase shift is 
7r (conventionally assumed for a bound state) and monotonically decreases 
to zero as the energy increases. Again we see a very important difference 
between the relativistic and nonrelativistic expressions for large (3\ values 
{Pi > m K ). At fix = 2000 MeV and E = 1400 MeV the difference is as large 
as 30°. Let us also notice that at higher Pi the phase shift decrease with 
energy is much less steep than at lower Pi = 150 MeV. This decrease at the 
KK threshold is governed by the negative value of the KK scattering length. 
The general expression for the KK scattering length in the relativistic case 
is given by Eq. ( |38|) provided we substitute in it the KK channel parameters 
mi, Pi and An. 

The values of the scattering lengths and the coupling constants are given 
in Table |. As indicated in the last column, the relativistic corrections at the 
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KK threshold amount to 13% if the value f3\ is as large as 2000 MeV, while 
at low f3\ < 150 MeV they are smaller than 1%. 

Next we examine the main difference between the antibound state and 
the bound state discussed above. The corresponding Jost function zero is at 
k\ = iai, but now a\ is negative. The nonrelativistic relation ([43"D leads to a 
weaker coupling than in the previous case (— 1 < An < 0). In Fig. lb we show 
the energy dependence of the phase shifts calculated in the relativistic and 
nonrelativistic cases for two values of the range parameters (3\ as in Fig. la. 
The phase shifts increase at the KK threshold starting from the zero value, 
which means that the scattering length is positive. Again we notice that the 
relativistic effects are mainly important for the large values of the parameter 
(3\ and for higher energies. The phase shifts calculated relativistically are 
larger that the corresponding nonrelativistic values. 

A special case in the KK channel is the existence of a resonance at the 
complex energy Ms — iTs/2, where Ms is the resonance energy and Ts is its 
width. In the complex momentum plane there are two S'-matrix poles in the 
lower half plane at k\ = ±/cr — ikj, where kj > 0. They coincide with zeros 
of the Jost function: 

D K (±k R - ikj) = 0. (44) 

For our nonrelativistic and simple choice of the KK interaction (Eq. ||) 
the interaction strength and the range parameter are fixed, since f3\ = kj and 
An = k\jk\ > (repulsive interaction). In the relativistic case, however, 
there are two possibilities of choosing the parameters f3\ and An at a given 
resonance position. One set of parameters is very close to the nonrelativistic 
set. This is related to the values of ki which in general must be smaller 
than about 100 MeV if we would try to attribute the discussed resonance 
to the observed rather narrow / (975) meson strongly coupled to the KK 
channel. The second relativistic solution is obtained for very large f3i value 
and An < (attractive interaction). The phase shifts are plotted in Fig. lc 
and their behaviour is very different in the two cases. For smaller (3i values 
close to 70 MeV (upper curve) we at first observe an energy decrease and then 
an increase of the 5 K -g function. This means that the resulting scattering 
lengths are negative. The phase shifts calculated for the relativistic and 
nonrelativistic propagators are very similar. They differ by less than one 
degree so the corresponding lines in Fig. lc are indistinguishable. The lower 
curve corresponds to the second relativistic solution and the phase shifts 



13 



increase staring from 0° at the threshold. In this case the scattering length 
is positive. 

At the end of this section we return to the discussion of the bound state, 
especially its wave function in the momentum and configuration spaces. An 
important parameter describing the bound KK system near the threshold 
is the root mean square radius parameter. It has been discussed in the 
nonrelativistic model [18], and its preferable value was given as (r^) 1 / 2 = 0.76 



fm. The nonrelativistic form of the radial wave function is simple and can 
be found in Ref. | lq| . Its Fourier transform ^n{p) describes the momentum 



distribution of the KK relative motion (see Ref. pl[|) 

^n(p) = c N G N (p) gi (p), (45) 

where Gn(p) = — ^k/^i +p 2 ) is the nonrelativistic propagator correspond- 
ing to the binding energy E R = a\/mx and cat is the normalization constant 

c N = -\J^{a 1+ ^ (46) 

such that J d 3 k | ^/at(/c) | 2 = 1. Similarly in the relativistic case we write 

* R (p) = c R G R (p) gi (p), (47) 



where G R (p) = l/(E — 2Jp 2 + m K ) is the relativistic propagator correspond- 



ing to the total energy E = 2y m 2 K — a 2 and c R is the normalization constant. 
The relativistic wave function in the configuration space 3> R (r) is related to 
ty R (p) by the Fourier transform. An important parameter describing the 
radial extension of the KK system is the root mean square diameter 

(r 2 ) = J d 3 rr 2 | $ fl (r) | 2 (48) 

or the root mean square radius (r 2 s ) l l 2 = ^(r 2 ) 1 / 2 . The results of numerical 
calculations are given in Table [Q]. The first set of parameters a\ and (3\ is 
taken from Ref. [fT8]] . The second and third line correspond to the results 
of fits described further in Sect. V. The root mean square radii are small, 
typically about 0.7 fm, so the KK molecule has a rather compact structure. 
This fact can have important consequences. For example, in Ref. |2T5j a value 
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1.2 fm has been used to predict the width of the radiative <fi — > fo (975)7 decay. 
If instead of 1.2 fm we use the value 0.7 fm the corresponding width increases 
by a factor of 2. Let us mention that the annihilation probability of the KK 
molecule into two photons will be substantially increased for small KK radii 
(compare |42]|). In Table II we can also notice a general trend: the root mean 
square radii calculated for the relativistic propagator are smaller by about 
0.1 fm in comparison with corresponding values for the nonrelativistic case. 
The difference comes from the behaviour of the radial wave functions for the 
values of r smaller than about 0.7 fm, where the relativistic wave function 
takes higher values than the nonrelativistic one. 



V 7T7T and KK coupled channel analysis 



In this section we describe the coupled channel analysis performed for the 
interacting irn and KK pairs. Starting from the irn threshold we put the 
upper energy limit at 1400 MeV, where one can still neglect the coupling to 
other channels with higher thresholds (as discussed in |19|). The interactions 
(ihH) i n both channels have the forms as simple as possible. As it was men- 
tioned in Sect. Ill the rank-two potential allows us to describe a substantial 
increase of the irn phase shifts below 1 GeV and the existence of the known 
/o(1400) meson. All the potential parameters have been obtained by fitting 
the calculated phase shifts 67m, 5 K j^ and the inelasticity parameter r\ to the 
experimental data. Above the KK threshold a sum ip = 6^ + 5 K j^ is often 
used and a quantity x = (1 — ?? 2 )/4 is introduced to represent the inelasticity. 
Except for the im threshold data of Refs. |p8| , [p9fl we have used the data of 

m 



Refs. pOl ) [pi]] covering a wide energy band. The total number of the Sqm 
points is 56. In addition, 17 points of the 77 dependence have been taken from 
Ref. [fUjfl . For the <p dependence we use two distinct data sets ||3 



44 which 



essentially differ in their behaviour near the KK threshold. In the data ffij 



used in the analysis of Ref. J44J a constant increase of ip is seen while in the 
data of Ref. some threshold decrease with energy can be observed. As it 
was shown in Sect. IV, this difference could be crucial for understanding the 
nature of /o(975), so we have performed two separate fits to the collection of 



the Smi, 5 K x and rj points. The data set containing 16 (^-points of Ref. [43 
is further called set 1 while in the set 2 we include 17 ^-numbers read from 
Ref. JPj] (the total number of the data points is 89 and 90 respectively for 
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these two data sets). 

The fitted parameters are shown in Table III. The asymmetric errors A + 
and A„ correspond to an increase of the total x 2 value by one unit. Let us 
notice that all the channel coupling constants An are negative, which means 
that the interactions are attractive in both channels. The couplings between 
channels are small. The ranges of the form factors are short since the /3j 
values are in the GeV region. During the fitting procedure we have found 
that it was convenient to use as an independent variable the product A 33 /3 3 
instead of /3 3 . For both data sets this product is very close to —7im n ~ —431 
MeV. This fact is not accidental as it is explained in Appendix. The very 
large value of /3 3 is compensated by very small values of the corresponding 
coupling constants A 33 and Ai 3 . The starting parameters A 2 2 and in 
the minimization procedure have been given in Sect. III. Other coupling 
constants were initially put equal to zero and $ parameters had the value 
1 GeV. In Table [TV] we show x 2 values of two fits for separated data parts. 
The data on Sjctt and rj are very well described by both fits but the data on 
(p only by the fit to the data set 1. 

In figures 2-4 the results of our fits to the data sets 1 and 2 are shown. 
Theoretical 5im and r] values differ only slightly for both fits. In Fig. 3 we 
notice an initial decrease of the ip curves at the KK threshold, so the data 
set 1 is clearly favoured by the model (compare also the corresponding x 2 
values of Table fV\) . Due to a very similar monotonic increase of the 6%% 
phase shifts above 1 GeV, the differences in the ip behaviour for the data sets 
1 and 2 are related to various possible trends of the S K ^. As discussed in 
Sect. IV the trend seen in the yj-data of Ref. JE| is connected with a presence 
of the quasibound KK state but in order to exclude other possibilities new 
precise S K -^ measurements near the threshold are needed. 

We have investigated a pole structure of the S matrix elements and the 
results are shown in Table [V]. A very wide / (500) meson (see Sect. Ill) 
causes a systematic increase of the 5^^ phase shifts starting from the tctt 
threshold. The presence of the narrow /o(975) manifests itself by a strong 
jump of the 5-jm and the enhancement of the inelasticity function x near 
the KK threshold (see Figs. 2 and 4). A further increase of the and tp 
above 1.2 GeV is related to the presence of the /o(1400) meson. Also the 
structure of the inelasticity around the energy 1.3 GeV is caused by this 
state. We should point out that the values of the resonance parameters may 
be changed if a fit to another data sets is performed. We expect that changes 
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may be much larger than the error ranges quoted in Table [V]. This is due to 
the fact that some experiments on the nir and KK scattering have supplied 
contradictory results, as was illustrated by our choice of the sets 1 and 2. 
We think, however, that the differences between data sets are representative. 
In particular the resonance parameters are not very much influenced if we 
compare results of two fits in Table |V|. Let us mention here that a moderate 
width of the / (1400) is obtained together with the presence of the wide 
/o(500) meson. The calculated masses and widths of the /o(1400) meson 
are in very good agreement with the corresponding values obtained in ref. 
[gH (M = 1425 ± 15 MeV and T = 160 ± 30 MeV). In the absence of the 
former state the systematic increase of the nir phase shifts below the energy 
of 1 GeV should be caused by a very wide / (1400) state (see for example 



Ref. |18|1 ) or by another wide state at the energy of about 1 GeV. 

In Fig. 5 the matrix element singularities are presented using the 
complex variable z (compare [IE] and [15] ) 



* = -4i^.. (49) 

The positions of poles and zeros have been calculated from the fit to the data 
set 1. A very similar structure of the S^nr has been found for the fit to the 
data set 2. The energy sheets on the z-plane are defined by the signs of the 
imaginary parts of the &2 and k\ momenta (Im A)2, Ith k\) as follows: I (+,+), 
II (-,+), HI ( - , - ), IV (+,-). In Fig. 5 we notice two poles on sheet II labelled 
as 1 and 2. They correspond to the /o(500) and /o(975) resonances. The 
/o(1400) resonance position on sheet III is indicated by 3. Those poles are 
the nearest singularities lying close to the physical region. The tctt threshold 
region is also strongly influenced by the cuts located on the imaginary z- 
axes. In the momentum space these cuts are on the fcj axes (from ± im,i to 
±ioo). Their origin is due to the presence of the logarithmic and square 
root functions in the Jij integrals (see Appendix). In the nonrelativistic case 
these integrals have much simpler structure and there are no cuts in the z— 
plane. Two poles at Kez > are particularly close to the | z |= 1 circle: 
pole 2 associated with the /o(975) resonance and the another pole lying on 
sheet III with a zero almost superposed on its top. Pole 2 can be related to 
the antibound state discussed in Sect. IV. If the coupling between channels 
is switched off then the poles move to the \z\— 1 circle. The positions of 
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other zeros or poles are connected with the discussed /q(500) and /o(1400) 
resonances. In particular the structure of the inelasticity function shown in 
Fig. 4 is a result of an interplay of the S'-matrix zeros and poles. The S-jm 
matrix element also has two second order zeros related to the form of the irir 
potential. They correspond to the values of k 2 = and k 2 = —i/3 3 and 
due to large values of /3 2 and /5 3 two of them lay very close to the origin of 
the z-plane and the other two very far from that point. These poles are not 
shown in Fig. 5. 

We have also performed fits with only one term in the 7T7T potential i.e. 
with A33 = in Eq. (f|). The x 2 fit was poor. In this case the SVtt matrix 
element has no pole corresponding to the / (1400) meson and the S^n and 
ip phase shifts have a very flat behaviour above the energy 1.2 GeV. 

Similar fits as for the relativistic form have been done using a nonrela- 
tivistic form of the Jost function. The coupling constants A 22 and A 33 had 
then opposite signs and the / (975) meson also appeared as a KK bound 
state. We have found that apart from the /o(975) meson there was only one 
very wide (r « 1300 MeV) state at the energy of about 850 MeV. The x 2 
value for the fit was 283 for 89 points so it was much worse than in the rela- 
tivistic case. Energy dependence of the 6%% , ip and x was similar to that one 



in the relativistic fit with only one term in the Tin potential. Like in Ref. QlS 
the S matrix also has two poles on the negative part of the imaginary axis 
k 2 but no poles corresponding to the / (500). 

Scattering lengths a^it and a K -^ are very important quantities describ- 
ing the near threshold tttc and KK interactions. Comparison of the values 
obtained by us with experimental measurements and other theoretical esti- 
mations is shown in Table [V|. Good agreement of the predicted a-jm values 
for the relativistic case with other data is found. In the nonrelativistic case 
there are two S'-matrix poles on the k 2 imaginary axis and one of them is 
relatively close to the tttt threshold, which gives a large positive contribu- 
tion to the scattering length. Therefore the scattering length is much larger 
than in the relativistic case. An appearance of the imaginary part of the 
KK scattering length is related to the fact that above the KK threshold the 
annihilation process into two pions is possible. The value of this imaginary 
part is evaluated using a low momentum approximation of the inelasticity 77 

77 ^l-2bk 1 + 0(k 2 1 ), (50) 



where the parameter b is equal to Ima^-^ (see also Ref. |nj). The negative 
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value of the real part of the KK scattering length is caused by an influence of 
a single pole lying close to the KK threshold on the sheet II (see Fig. 5). This 
pole is related with the narrow / (975) meson. We expect that in a case when 
two poles are close to the KK threshold on sheets II and III (as discussed 
in Ref. ||17|| ) the absolute value of the real part of the KK scattering length 
can be much smaller. The reason is that the contributions of such poles have 
opposite signs so they cancel to a large extend. The cancellation is complete 
if the poles lie symmetrically in the kaon momentum space (at k\ and -k\). 
The precise experimental measurements of this quantity could distinguish 
between these two possibilities. There is one experimental estimation of the 
KK scattering length based on the data of ref. |47j which gives numbers 



similar to our findings. Unfortunately, the KK phase shifts used by Wetzel 
et al. do not satisfy a very general requirement S K j^(— ki) = —S K j^(ki), so 
the values of the scattering lengths obtained in [M could be questioned. 
The K + K~ scattering lengths quoted in ref. [IS] and given in Table [VT] are 



theoretical estimations used in the calculations of the properties specific to 
the kaonic atoms. Contrary to the KK case, more information on the tttt 
scattering length is available. Our values for the relativistic fits are in good 
agreement with other determinations (see Table |VID . 

Due to the channel coupling the /o(975) and /o(1400) mesons can decay 
into the tttt and KK pairs. The / (975) meson being presumably a KK 
bound state becomes unstable due to the annihilation process KK — > tttt. 
The branching ratio of its decay into the KK channel is 

B = ^, (51) 

J- tot 

where T tot denotes the total and T K j^ the KK fractional decay width. If we 
wish to compare the experimental /o(975) branching ratio with the theoreti- 
cal one, we must take into account the fact that due to the proximity of the 
KK threshold the /o(975) resonance has no a typical Breit-Wigner shape. 
As it was discussed in Ref. |TS| the average branching ratio B av of the / (975) 
reads 

j C M m ax 

B ™ = Wm m L dE{l - v ) ' (52) 

l{M max - M s ) JM th 

where M max , M s and M t h denote respectively the upper limit of integration, 
mass of the meson and mass of the KK threshold. In the numerical cal- 
culations we put M max = 4T to t- We obtain B av = 18.0 ± 0.7% for the fit 
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to the data set 1 and B av = 18.4 ± 0.6% for the fit to the data set 2 in 



agreement with the Particle Data Group value 21.9 ± 2.4% j34j. In the case 
of the /o(1400) meson we have performed the integration from M s — T tot to 
M s + Y tot substituting 2Y tot in place of 2(M max — M s ) in Eq. fl5"2"|). For the 
fo(1400)^KK branching ratio we obtain 16±1% for both data sets in agree- 



ment with [43]. Let us notice that the /o(1400) resonance mass is slightly 



higher than the upper energy limit 1.4 GeV used in the fits. This fact may 
lead to some uncertainties in the predicted values of the /o(1400) branching 
ratios. The Particle Data Group value is 6.4^5% @- ^ ^ s entirely based 
on the analysis done in Ref. in the energy range from 1100 to 1420 MeV. 
The corresponding value of the branching ratio calculated for the parameter 
set 1 in the same energy range is as high as 30.3%. This is a result of the 
existing discrepancy between the s-wave intensities obtained by the Argonne 
group E5| and the CERN- Cracow-Munich group [PJ. In our analysis we 
have chosen the Argonne data, since they include points closer to the KK 
threshold. For energies higher than 1.4 GeV an influence of other scalar- 
isoscalar states such as /o(1525), /o(1590) or even / (1710) (see references in 



5j]) may be important. We have also tried to include in our analysis the 
energy range from 1.4 GeV to 1.5 GeV and we could not simultaneously fit 
the Sun and if phase shifts. The 5jctt data in this energy region favour wider 
/o(1400) meson; the sum if = 5^^ + S K -^ which increases steeply with energy 
requires an existence of a narrow state. 



VI Conclusions 

We have performed the analysis of the isoscalar spin zero tttt and KK 
coupled channel interactions. The model was based on a separable potential 
formalism. Two channel scattering amplitudes have been evaluated from the 
coupled equations of the Lippmann-Schwinger type. Using the relativistic 
propagators we have obtained very good fits to the experimental tttt and 
KK scattering data. It has been found that the potentials in the inr and 
KK channels are attractive and have short range form factors. The y 2 fits 
in the relativistic case were 3 or 4 times better than in the nonrelativistic 



case. In contrast to the nonrelativistic calculations (see QT8[ and Sect. V) new 
/o(500) state has been found in the irn interactions treated relativistically. 
Its large width about 500 MeV is in agreement with other estimations. In the 
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relativistic version of our model the mass and width of the /o(1400) meson 
are well described (see Table 0). Relativistic effects also play the important 
role near the 7T7T threshold. The tttt scattering length is significantly smaller 
than in the nonrelativistic case and is in agreement with the experimental 
data and other models (see Table [V| ). 

We have also been interested in describing the KK interaction. Two 
different data sets on the KK phase shifts have been taken into account. 
Our solutions indicate the existence of a quasibound state below the KK 
threshold which we identify with the /o(975) meson. Thus this state does 
not seem to be a typical qq meson. Some observables, e.g. the KK scattering 
length and root mean square radius of the /o(975) have been evaluated. We 
predict the large negative value of the KK scattering length about —1.7 fm 
and the corresponding imaginary part of the order of 0.6 fm (Table VI). In 
our calculations the /o(975) state appears as the compact KK system with 
the root mean square radius of about 0.7 fm, so it is not a KK molecule of 
the deuteron size. Analysis of the KK channel leads to a conclusion that at 
the KK threshold the relativistic corrections to the KK scattering length 
or the root mean square radius are of the order of 10% and they gradually 
grow with increasing energy. 

The existing data above the KK threshold are still controversial. Com- 
parison of our expectations with new experimental data may be done in near 
future and a new light on the nature of the /o(975) and ao(980) mesons can 
be shed on if very precise measurements near the KK production threshold 
like those planned at COSY [|5(J] are performed. Very good energy resolution 
would enable us to take into account a mass splitting between the charged 
K~K + and neutral K°K° modes. Due to different interactions of the K~ 
and K + with nucleons in the GeV region the molecular picture of the /o(975) 
can also be verified by studying its interaction inside nuclear matter (see also 



Ref. p|). 
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Appendix 



In this appendix we give explicit expressions for the t-matrix elements 
( |TTD in terms of the coupling constants A^ and the integrals iy (Eqs. (UP 
and (0)): 



'11 

tl2 
tl3 
t 2 2 

^23 

^33 



+ 



D 1 [An — (A11A22 — A^ 2 )/22 — (A11A33 — A^ 3 )/33 + 2A12A13/: 

4)], 
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(A11A22A33 — A33A^ 2 — A22Ai 3 )(/22-^33 
= D X (Ai2 — A12A33/33 + A22Ai3/23); 
= D (A13 — A13A22/22 + A33A12/23), 

— D 1 {Ai2/n[Ai2(l — A33/33) + A22A13/23] 

+ A22(l — Aii/n)(l — A33/33) — A22Ai3/n(Ai2/23 + A13/33)}, 

— D 1 [A22A 3 3/23 + A12A13/11 

— (A11A22A33 — A22A13 — A33Ai 2 )Al-^23] 5 
= D 1 {Ai3/n[Ai3(l — \22h2) + A33A12/23] 

+ A33(l — Aii/n)(l — \22h2) — A33Ai2/ll(Ai3/23 + A^-^)}- 



(Al) 
(A2) 
(A3) 

(A4) 

(A5) 

(A6) 



In Eqs. ( |A1[ )-(|A6|) D denotes the function defined in Eq. (|24|). The integrals 
( [14]) and ( [15]) have also been used in |23| to calculate the K~-N poten- 
tial but analytical formulas have been given only when the propagator G in 
Eq. (|2]) was nonrelativistic (see also |)lj). We show below that in the rel- 
ativistic case the analytical calculations can be performed as well. A very 
helpful substitution is: y = y 1 s 2 /(s 2 + mf). Then after some straightforward 
algebraic operations the integrals ly- can be rewritten in the terms of the 
integral 



Fix 2 



dy 



o 1 — y 2 x 2 + ie 



(A7) 



If x is real then 



F{x 2 ) 



1 

Ix 



171 

"Zx~ 



arctan(| x \)J \ x 



if x 2 > 1, 
if < x 2 < 1, 
if x 2 = 0, 
if x 2 < 0. 



(A8) 
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Finally the integrals can be expressed in terms of Jjj (Eq. (|T9|) ) : 



E; 



4m,- fl 



7^(1 + a?) 2 



■i(l + «?)^ + a* a C?' 



(A9) 



^23(^2 



+ 



2m 2 (l - za 2 )(l - z"a 3 ) ft + A Trv^M^l + a|)(l + a 



4ftft 2 ft 2 -ft 2 



rf 3 (l + a 2 )(6 2 -l)-rf 2 (l + a 2 )(6 2 -l) . (A10) 



In (P) and ([ATo| ) 



E; 



2Jkf + ml a 



Ei_ 
2h 



Hi = F{Cl 



and 



1,2,3. 



(All) 



(A12) 



In ( |A11|) and ( |A12j ) the m,i are defined by (0) and k 3 = k 2 . 

The nonrelativistic limit of the integrals ( |T4] ) and ( |T5D can be obtained 
from (|A9|) and ( |A10| ) when m ; — > 00. Physically it means that the nonrela- 
tivistic limit is achieved when both hi [i = 1, 2) and ft (j = 1, 2, 3) are much 
smaller than m;. Then 

= — TTv? (A13) 



and 



^23(^2 



2 v / ftft 



(A14) 



(/?2 + /?s)(l - «a 2 )(l - ia 3 ) 

Let us discuss the limit ft — > 00. For very large (3% the integral J33 is 
proportional to /5s: 

ft 



J: 



33 



vrm,, 



(A15) 



In order to compensate too large values of J33 in Eq. (26-^8) we can multiply 
it by a sufficiently small A 33 value which should be inversely proportional to 
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/3 3 . If 03 — > oo and A 33( 5 3 = —nm^ then the term 1 — A 33 J 33 in Eq. (f26|) 
does not grow with /3 3 . 

In the limit /3 3 — > oo the integral J23 vanishes in Eq. (^) and 
the Jost function Dti-(/c 2 ) is approximately a product of two terms 
(1 — A22J22XI — A 33 J 33 ). The zero of the first term can be attributed to 
the /o(500) meson as discussed in Sect. III. For the values /3 3 and A 33 as 
given in Table III the second term (together with the small Jf 3 term) allows 
us to obtain the second zero relatively close to the real axis /c2- In Table [V] 
we assign this zero to the /o(1400) meson. 
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Table I: KK coupling constants An and scattering lengths a corresponding 
to the KK bound state at 975 MeV (without coupling to the nn channel). 



A 


nonrelativistic 


relativistic 


1 M. 

1 a N 




(MeV) 


An 


a N (fm) 


A n 


a R (fm) 


- 1 | -100% 


150 


-2.55 


-4.32 


-2.48 


-4.28 




1 


500 


-1.39 


-2.81 


-1.18 


-2.70 




4 


2000 


-1.09 


-2.35 


-0.54 


-2.05 
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Table II: Root mean square radii of the KK wave function. 



E 
MeV 


a\ 
MeV 


A 

MeV 


nonrelativistic 

(r|)V2 

fm 


relativistic 

(r|)V2 

fm 


971.65 


98.90 


2000 


0.76 


0.66 


973.36 


94.51 


1496 


0.81 


0.73 


973.71 


93.59 


2177 


0.80 


0.69 
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Table III: Model parameters and their up (A + ) and down (A_) errors. 



fitted 


set 1 


set 2 


parameters 


values 


A + 


A_ 


values 


A + 


A_ 


An 


-0.658 


0.030 


-0.031 


-0.511 


0.075 


-0.078 


A22 


-0.201 


0.003 


-0.003 


-0.201 


0.004 


-0.004 


A 33 x 10 5 


-7.46 


2.32 


-3.51 


-8.95 


3.11 


-5.55 


A12 


0.0363 


0.0021 


-0.0025 


0.0251 


0.0064 


-0.0048 


Ai 3 x 10 6 


3.0 


2.5 


-1.0 


2.8 


2.9 


-0.9 


A (GeV) 


1.496 


0.115 


-0.082 


2.177 


0.528 


-0.402 


fa (GeV) 


1.162 


0.052 


-0.051 


1.141 


0.052 


-0.052 


(3 3 xA 33 (MeV) 


-431.162 


0.070 


-0.040 


-431.145 


0.065 


-0.035 


(3 3 x 1(T 3 (GeV) 


5.8 


2.6 


-1.9 


4.8 


2.6 


-1.8 



Table IV: \ 2 values for the fits to the data sets 1 and 2. 



Set 


x 2 


x 2 


x 2 


x 2 


No 


total 


dim data 


if data 


rj data 


1 


75.2 


52.2 


6.7 


16.4 


2 


100.8 


50.5 


34.2 


16.1 
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Table V: Masses and widths of resonances obtained in fits to the data sets 
1 and 2 compared with values of the Particle Data Group [34] and Ref. [39] 
for / (500). 



pole 


set 1 


set 2 


Particle Data Group 


M (MeV) 


T (MeV) 


M (MeV) 


T (MeV) 


M (MeV) 


T (MeV) 


/o(500) 
/o(975) 
/o(1400) 


506±10 
973±2 
1430±5 


494±5 
29±2 
145±25 


505±10 
974±2 
1428l^ 3 


497±5 
30±1 
15711 


< 700 
974.1±2.5 
~ 1400 


> 600 
47±9 
150 400 



Table VI: Comparison of the tot and KK scattering lengths obtained in the 
present work (sets 1, 2 and nonrelativistic fit) with other determinations. 



Source 




Re(a x¥ ) (fm) 


Im(a^) (fm) 


set 1 


0.172±0.008 


-1.73±0.07 


0.59±0.04 


set 2 


0.174±0.008 


-1.58±0.09 


0.61±0.04 


nonrel 


fit 


0.40 


-1.70 


0.78 




48 






-1.15 


1.80 




32 




0.207±0.028 
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FIGURE CAPTIONS 



Fig. la. KK scattering phase shifts versus energy in the presence of the 
bound state at 975 MeV. Solid lines correspond to the relativistic KK 
propagator, while the dashed line to the nonrelativistic propagator. 
Two upper lines are calculated for /3i=2000 MeV, two lower lines - for 
#l=150 MeV. 

Fig. lb Same as in Fig. la but for the antibound state at 975 MeV. 

Fig. lc KK scattering phase shifts versus energy in the presence of the 
resonance at M s = 993 MeV of the width T = 46 MeV. Both curves are 
calculated for the relativistic propagator: the upper one corresponds 
to p!=72 MeV and An=1.16 and the lower one to /3i=10.47 GeV and 
An=-0.14. 

Fig. 2. Isospin s - wave tttc phase shifts calculated for the data set 1 
(solid line) and set 2 (dashed line). Arrows indicate the energy range 
used in the data fits. Data are from [28-31]. 

Fig. 3. Sum of the irir and KK phase shifts. Full circles denote data 
set 1 [H3|, open circles - data set 2 [04]]. Lines and arrows as in Fig. 2. 



Fig. 4. Inelasticity parameter x = (1 — 7/ )/4. Data are taken from [j43 
Lines and arrows as in Fig. 2. 

Fig. 5. Structure of the SVtt matrix element in the complex z - plane. 
Positions of poles (1 - / (500), 2 - / (975), 3 - / (1400)) are indicated 
by crosses and zeros by circles. The roman numbers label the energy 
sheets. The bold line shows the physical region and the rectangles on 
the imaginary axis indicate the S^k cuts. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403264vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403264vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403264vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403264vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403264vl 



